function psiout = draw4triangle(Nx,Ntop,psiin)
%DRAW4TRIANGLE draw a ST with scale L=4 & top # Ntop
if mod(floor(Ntop/Nx),2)==1
    [psiin(Ntop),psiin(Ntop+Nx),...
     psiin(Ntop+Nx+1),psiin(Ntop+2*Nx-1),...
     psiin(Ntop+2*Nx+1),psiin(Ntop+3*Nx-1),...
     psiin(Ntop+3*Nx),psiin(Ntop+3*Nx+1),...
     psiin(Ntop+3*Nx+2)]=deal(1);
else
    [psiin(Ntop),psiin(Ntop+Nx-1),...
     psiin(Ntop+Nx),psiin(Ntop+2*Nx-1),...
     psiin(Ntop+2*Nx+1),psiin(Ntop+3*Nx-2),...
     psiin(Ntop+3*Nx-1),psiin(Ntop+3*Nx),...
     psiin(Ntop+3*Nx+1)]=deal(1);    
end
psiout=psiin;
end

